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The flavor singlet axial and tensor charges of the nucleon are calculated in lattice QCD. 
We find AE = Am + Arf + As = +0.638(54) - 0.347(46) - 0.109(30) = +0.18(10) for the 
axial charge and SJ: = Su + Sd + Ss = +0.839(60) - 0.231(55) - 0.046(34) = +0.562(88) 
for the tensor charge. The result for the axial charge shows reasonable agreement with 
the experiment and that for the tensor charge is the first prediction from lattice QCD 
' before experimental measurements. 

^ ' 1. Introduction 

The parton spin structure of the nucleon at the twist two level is characterized by two 
^ ■ structure functions gi{x, /i) and hi{x, fi) with x being the Bjorken variable and /i being the 
renormalization scale (see Ref. The functions gi and hi represent the quark- helicity 
^ ! distribution and the quark-transversity distribution respectively. Experimentally the for- 
^ [ mer has been measured by deep inelastic lepton-hadron scattering(DIS) experiments since 
1980 's. On the other hand, we have no experimental data for the latter, which could be 
measured in polarized Drell-Yan processes as well as in inclusive pion or lambda produc- 
tions in DIS. The polarized Drell-Yan experiment is planned using the relativistic heavy 
ion collider (RHIC) at BNL||. 

The first moment of gi is called the axial charge of the nucleon Aq, which is related 
to the nucleon matrix element of the axial vector current through the operator product 
expansion. The flavor singlet component of the axial charge AS = Au + Ad + As, being 
interpreted as the total fraction of proton spin carried by quarks in the parton picture, 
has been a matter of great interest since results of the European Muon Collaboration 
experiment indicated that AE is quite small AS = 0.12(17) and that the strange quark 
contribution is unexpectedly large As = —0.19(6)0. Further experiments have since 
been performed with protonQ-^, deuteron[^-|^ and neutron |1^,|11]] targets. Recent global 



analysis of these experimental data gives AS = 0.10io}[ at the renormalization scale 
/i2 = oolTg. 

The first moment of hi is called the tensor charge of the nucleon Sq, which is written in 
terms of the nucleon matrix element of the tensor operator. It is of great importance to 
predict the value of 6q before experimental measurements. We are interested in whether 
the flavor singlet tensor charge 61^ = Su + 6d + 6s has also a small value as the flavor 
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singlet axial charge. 

In this report we first give a brief description of the calculational method for nucleon 
matrix elements in lattice QCD. We present our results for the axial and tensor charges 



including both the connected and the disconnected contributions [jT3|,|T^. For the axial 
charge our main concern is a direct check of the small value of AE and the large negative 
contribution of As. The result for the tensor charge, which is the first one from lattice 
QCD, is compared to that for the axial charge. The current status of the lattice QCD 
simulations for the flavor singlet nucleon matrix elements is summarized in Ref. |jl5|. 

2. Formulation of the calculational method of nucleon matrix elements 

The axial charge Ag and the tensor charge 5q are related to the nucleon matrix elements 
of the axial vector current and the tensor operator respectively, 

(P, ■sl^ZT/.Tsgb, s) = 2Ms^Ag, (1) 
(p, s|g2(7^^75g|p,s) = 2{s^py - s^p^) 5q, (2) 

where is the nucleon four momentum, M the nucleon rest mass and the nucleon 
covariant spin vector normalized as = —1. We will explain how to calculate these 
nucleon matrix elements using lattice QCD. Hereafter for the definiteness we consider the 
proton for the nucleon and concentrate on the case of zero spatial momentum. 

We consider a four- dimensional Euclidean space-time lattice with lattice spacing a, 
where sites are labeled hy x = {riia, n2a, risa, n^a) with ni, ■ ■ ■ , n4 integers. The interpo- 
lating field for the proton at site x = {x, t) is constructed to have the correct quantum 
numbers of the proton, 

Op^{x) = e"^" [{^u\x)C-^-i^d\x)) u\x) - {^(r{x)C-^-i^u\x)) u'{x)\ , (3) 

where s denotes the proton spinor, a, 6, c represent the color indices and C = 7472 is the 
charge conjugation matrix. The proton mass mp can be extracted from the exponential 
decay of the two-point function projected onto the zero spatial momentum state, 

C,{t) = (0| VOp,(f,t)Op,(0,0)|0) ^5^* (^\Ops\p,s){p,s\OpM ^_^^, ^ ^ ^ 
^ 2mp 

Unwanted contributions of the excited states fall off rapidly as t increases, because their 
masses are heavier than that of the proton. 

The local quark bilinear operator on the lattice is defined as 

0'^''ix)^qix)Tqix) iq = u,d,s) (5) 

where T = 7^75 (/i = 1, ■ ■ ■ , 4) is employed for the axial charge and T = a^y = [7^7^ — 
7y7^]/2 (/i, z/ = 1, ■ ■ ■ , 4) for the tensor charge. The three-point function projected onto 
the zero spatial momentum state is used to find the forward proton matrix element of the 
operator O^**; 

Csit) = (0|X:Oft(x,t)Eor(^',OOft(0,0)|0) (6) 

X x' 

2mp 2mp 
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We define the ratio of tlie tliree-point function to the two-point function, whose time 
dependence is expected to be 

{o\EnOpMt)OpMom 

large ^ , (p, S | Op " |p, s) . ^ ^, „ 

— > const. + ^ — -^—^ t ioi t> t > 0, 
2m p 

where the sum over t' is employed to increases statistics. Avoiding the excited state 
contaminations in small t region we can extract the matrix element (p, s|Op**|p, s) by a 
linear fit of i?Or (^) terms of t. 

We should note that the operator Op** defined on the lattice is different from the 
operator O^'^^ defined in the continuum regularization scheme, e.g., MS. The operator 
defined in each regularization scheme can be related by a renormalization factor Z^p, 

Or'^(/i) = Zo,(/ia, a,)Oi?"(l/a), (8) 

where /x denotes the renormalization scale and a is the lattice spacing. Using this renor- 
malization factor we can convert the nucleon matrix element on the lattice into that in 
the continuum MS scheme, 

(p,s|Or'^(/i)|p,s) = Zo,(/ia,a,)(p,s|Oi?"(l/a)|p,s). (9) 

can be estimated either perturbatively or non-perturbatively. In this report we use 
the perturbative estimate for Zq^ ||16| . 



3. Monte Carlo techniques for calculation of two- and three-point functions 

In the path-integral formulation the proton two-point function of is represented by 

C2{t) = ^ / n ^^^M(2/)^5»^9(^)Eop.(x,t)Op.(o,o)e-^^^^"^-s.'^^'^^''^^ (lo) 

where Sg{Afj_} is the gauge action on the lattice and qDq{A^}q denote the quark action 
for each flavor q. Time is rotated into the imaginary axis and Z is the normalization 
factor so that (0|1|0) = 1. We can evaluate this functional integral for A^, q and q using 
Monte Carlo techniques. However, fermions are represented by anti-commuting Grass- 
mann numbers in the path integral, which we cannot manipulate directly with numerical 
computations. Therefore we perform the integration for the fermion fields analytically, 

C2{t) = \ I n dA,{y) Y: Op^ (f , t)Op^ (0, 0) n det^J A4e-^«{^''> (11) 

^ y X q 

with the overline meaning the quark field contractions 



OpXx, t)OpM 0) = " ^ C: (x; 0)Tr{C-S5C (x; 0)^-^5*6^ (x; 0)} + ■ ■ ■ , (12) 

X 

where Gq{x] 0) denotes the quark propagator, which is defined as the inverse of the quark 
matrix D^jA^}, connecting the origin and a site x = {x,t). Since it takes prohibitively 
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(a) connected (b) disconnected 

Figure 1. (a) Connected and (b) disconnected contributions to the three-point function 
of the proton and the operator Or- 

large computer time to evaluate the fermion determinant detDg in ( p^O]) we neglect the 
contribution of this term putting detDg = 1 for all flavor q. This is the "quenched 
approximation" often employed in current lattice QCD simulations. In this approximation 
virtual quark loop effects are not taken into account. 

A brief sketch of the calculation of C2{t) with the Monte Carlo techniques is as follows. 



an ensemble of gauge configurations A^S' 



.Ng) is generated with the distribu- 



tion proportional to exp {—Sg{A^}). 

ii) solve the quark matrix inversion on the i-th configuration numerically, 



Y.Dg{A<^^}{y;x)G'^\x-0) = S, 



y,o- 



(13) 



iii) in terms of the quark propagator G^g \x;0) the i-th proton propagator projected to 
zero spatial momentum [J2ssOp^{x,t)OpX0,0)Y^^ is constructed. 

iv) an average over the ensemble gives the vacuum expectation value, 



C2(t) = lim — y 



N, 



9 i=l I 



n(0 



(14) 



The three-point function C3 of (^ has two types of quark field contractions, represented 
by a connected [Fig.|l|(a)] and a disconnected [Fig.[l|(b)] diagrams. We should note that 
in lattice QCD simulations one can treat both the OZI preserving (connected) and OZI 
violating (disconnected) contributions even in the quenched approximation. The calcula- 
tion of the connected diagram requires the quark propagator inserted with the operator 
Qiatt pj^ojected to zero spatial momentum at the fixed time slice t'] 



G«(x,t;0,0) 



J2G^^\x,t;x',t')TG^^\x',t';6,0). 



(15) 



We can obtain this quark propagator by solving |T7 
ED,{4^}{y;x)G^;\x;0) = ^5,,.,rG»(x';0). 



(16) 
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There exists a serious difficulty in the calculation of the disconnected piece . To project 
the operator with zero spatial momentum we need the sum of quark loops J2x' t'; x', t') 
at the time slice t'. The conventional point source method of (0) for each source {x',t') 
requires quark matrix inversions, which would take about a few months of computer 
time for each configuration if we employ standard parameters used in current numerical 
simulations. To overcome this difficulty we use a variant of the method of wall sources |]TK|, 



which makes use of local gauge invariance of QCD. We prepare a quark propagator solved 
with unit source at every space site at the time slice t' without gauge fixing: 

= (17) 

where G^^\x,t) = Y.x' G''g\x,t; x' ,t') from the linearity of this equation. The product of 
the nucleon propagator and J2x'~^A^G''g\x,t')] equals the disconnected amplitude up to 
gauge variant non-local terms which cancel out in the average over gauge configurations 
due to local gauge invariance. The superior feature of this method is that it requires 
only one quark matrix inversion, instead of inversions for each gauge configuration to 
calculate the spatial sum of quark loops. For the three-point function in the numerator 
of Ro^{t) of (^, where we take the sum of the operator over t', extension of the above 
calculational methods for the connected and disconnected amplitudes is straightforward. 



4. Parameters of numerical simulation 

For the lattice action in ( PH]) we use a single plaquette gauge action and the Wilson 
quark action whose parameters are P = G/g"^, g being the coupling constant, and the 
hopping parameter K. The parameter j3 controls the lattice spacing a and the hopping 
parameter gives the lattice bare quark mass through m^a = {l/Kg — for each 

flavor g, where K^. is the critical hopping parameter at which the measured pion mass 
vanishes. 

In Table |l| we summarize our simulation parameters. Our calculation is carried out in 
quenched QCD at /3 = 5.7 on a 16^ x 20 lattice. The up and down quarks are assumed 
to be degenerate with = = K^d, while the strange quark, which emerges as the 
quark loop in the disconnected diagram, is assigned a different hopping parameter Kg. 
We employ three hopping parameters i^^^, i^^s=0.160, 0.164 and 0.1665 corresponding to 
the physical pion mass ~ 970, 720 and 520MeV respectively, where the lattice spacing 
a is estimated from m^a in the chiral limit using rrip = 770MeV. The physical spatial size 
La is approximately 2.2fm. We analyzed 260 gauge configurations for the axial charge 
and 1053 configurations for the tensor charge, which are generated by the pseudo heat 
bath sweeps and are separated by 1000 sweeps. In order to avoid contamination from 
the negative-parity partner of the proton propagating backward in time, we employ the 
Dirichlet boundary condition in the temporal direction for quark propagators. We fix 
gauge configurations on the t = time slice, where the proton sources are placed, to the 
Coulomb gauge to enhance proton signals. Statistical errors for all measured quantities 
are estimated by the single elimination jackknife procedure. 

The strong coupling constant in the MS scheme is evaluated from the tadpole-improved 
one-loop relation between the MS coupling and the lattice bare coupling, aj^(7r/a)~^ = 
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Table 1 

Simulation parameters for the calculation of (a) the axial charge and (b) the tensor 
charge. rriNa is the nucleon mass in the chiral limit. 

13 = 5.7, X T = 16^ X 20, K = 0.16, 0.164, 0.1665 

^conf. Kc a^^{GeV) mjqa ay^{l/a) rriuda rUstt 

~(a) 260 0.16941(5) 1.458(18) 0.773(13) 0.2207 0.00340(9) 0.0826(22) 

(b) 1053 0.16925(3) 1.418(9) 0.8024(65) 0.2158 0.00347(5) 0.0896(12) 



Pav^iatt + 0.30928 with Pav the Monte Carlo expectation value of the plaquette[P^]. Re- 
ducing the scale from vr/a to 1/a by two-loop running we obtain a^{l/a), which is used 
for the calculation of the perturbative renormalization factors. An estimate of physical 
nucleon matrix elements requires the physical value of the degenerate up and down quark 
mass m„a = m^a = m^^a and that of the strange quark mass m^a. From fits of the 
measured meson mass spectrum we find (m^ra)^ = BpsmudO' and nipa = Ay -|- BymudCL- 
Using the physical vr to p mass ratio mT^/nip = 0.18 we obtain the degenerate up and 
down quark mass niudO'- The strange quark mass is estimated by generalizing the relation 
{m^aY = BpsTTiudO' to {iriKCiY — Bps{fnudO' + ^so) /"^ and using the ratio mK/mp = 0.64. 



5. Results 

5.1. Axial charges 

We calculate the ratio of (^ employing the axial vector current Ag = q'jp'y^q {q = 
u,d,s). In Fig. ^(a) we plot the connected contribution of u and d quarks to the ratio 
RAg{t) for the case of Kud = 0.164. A clear linear behavior is observed up to t ^ 14 for 
both contributions. For the disconnected contribution the quality of our data is not very 
good, in spite of high statistics of the simulation (see Fig. |^(b)). The region showing a 
linear dependence is limited to t ~ 5 — 10, and errors grow rapidly with increasing t; for 
t > 12 the signals are lost in the large noise. Nevertheless, we can still observe that the 
slope of the ratio for the disconnected contribution is smaller than that for the connected 
contribution and it has a negative value. To extract the axial charge for the connected 
and disconnected contributions we fit the data for Ra^ (t) to the linear form (^ with the 
fitting range chosen to be 5 < t < 10. 

The axial charges defined on the lattice are converted to those in the continuum MS 



scheme using the tadpole- improved renormalization factor[16,20 



^A,(aMs(Va)) = - [1 - 0.31aMs(V«)] • (18) 

The results for Am, Ad and As in the MS scheme are summarized in Table 0. We should 
note that the flavor singlet axial vector current requires an additional lattice-to-continuum 
divergent renormalization from diagrams containing the triangle anomaly diagram. We 
leave out this factor, since the explicit form of this contribution which starts at two-loop 
order has not been computed yet. 

We present the axial charges for u, d and s quarks in Fig. |^ as a function of the bare 
quark mass rriud in physical units. As we already remarked the values of the disconnected 
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Figure 2. Connected and disconnected con- Figure 3. Axial charges for u and d quarks 

tributions to R^^ (t) for u and d quarks at as a function of degenerate u and d quark 

Kud = 0.164. Solid lines are linear fits over mass. Open symbols denote the values ex- 

5 < t < 10. trapolated linearly to the chiral limit. 



contribution(circles) are small and negative. Their magnitude increases slightly toward 
the chiral limit, while the connected contributions decrease. 

We calculate the physical values of matrix elements in the following way. For u and d 
quarks, we estimate the sum of disconnected and connected contributions by first com- 
bining the two contributions in the ratio RAq{t) and then fitting the result to the linear 
form (|^ over 5 < t < 10 for each Ku^. The fitted values are extrapolated linearly to the 
chiral limit rriud = 0, where we neglect the degenerate u and d quark mass (see Table 0). 
For the strange quark contribution similar extrapolations to rriud = are made for each 
Kg, and their results in turn are interpolated to the strange quark mass m^. This analysis 
yields for the quark contribution to proton spin, 

AS = Am + Ac/ + As = +0.638(54) - 0.347(46) - 0.109(30) = +0.18(10). (19) 
These values, notably the sign and magnitude of the strange quark contribution, show a 



Table 2 

Proton axial charges in the MS scheme as a function of Kud- Addisc equals Audisc- Values 
at Kc are obtained by a linear fit in 1/K. 
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AUdisc 
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/s:,=0.1665 


0.1640 


0.1600 


0.1600 


0.9071(92) 


-0.2470(35) 


-0.025(10) 


-0.025(10) 


-0.035(12) 


-0.042(15) 


0.1640 


0.839(19) 


-0.2382(87) 


-0.066(23) 


-0.049(19) 


-0.066(23) 


-0.076(27) 


0.1665 


0.818(39) 


-0.231(23) 


-0.093(54) 


-0.068(35) 


-0.084(46) 


-0.093(54) 




0.763(35) 


-0.226(17) 


-0.119(44) 
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reasonable agreement with the phenomenological estimate AE = 0.10]'^o;iJ|jT2 . 

Possible sources of systematic errors in our results are scaling violation effects due 
to a fairly large lattice spacing a ~ 0.14fm aX (3 = 5.7, quenched approximation and 



uncertainties in the perturbative estimate of the renormalization factor (|T8|). For the first 
two systematic errors we can roughly estimate its magnitude from our result of the flavor 
non-singlet axial charge qa = Am — Arf = 0.985(25) at rriud = 0, which is about 25% 
smaller than the experimental value qa = 1.2601 (25) [^. The small value of qa-, possibly 
arising from these uncertainties, suggest that our result for AS might be underestimating 
the continuum value by a similar magnitude. For the perturbative renormalization factor 
we note that the lack of two-loop calculation for the flavor singlet lattice-to-continuum 
renormalization factor makes it difficult to specify the scale at which AS is evaluated, 
although we expect the scale dependence to be weak, being a two- loop effect. These 
points should be examined in future studies. 

5.2. Tensor charges 

So far, there exist several model calculations of the tensor charge. Non-relativistic 
quark model predicts 6u = Au = 4/3 and 6d = Ad = —1/3, while relativistic quark 
wave functions with non- vanishing lower components lead to 6u = 1.17 and Sd = —0.29 
together with the inequality \Sq\ > \Aq\ (g = m, (i)[^. There also exist attempts to 
calculate 6q using QCD sum rules {Su = 1.33(53), 6d = 0.04(2)) P^ , pB| and a chiral quark 



model {Su = 1.12, 6d = —0.42,5s = — 0.008)|Q. The main deficiencies of these model 
calculations are as follows. 

i) the renormalization scale where the matrix elements are evaluated is not clear. The 
tensor current has an anomalous dimension at the one-loop level, while that of the flavor 
singlet axial current starts from the two-loop level. 

ii) it is hard to estimate the contribution of OZI violating process, especially strange quark 
contributions, in a reliable manner. 

Lattice QCD calculation is free from these problems even in the quenched approximation. 

We calculate the ratio of (|^ for the tensor operator Tg = qa^uQ (g = u, d, s). The ratio 
RTqit) are plotted in Fig. || for the case of Kud = 0.164. Fig. ^(a) shows a clear linear 
behavior for the connected contributions up to t ~ 13, beyond which errors grow rapidly. 
For the disconnected contribution in Fig. H(b) the data stay around zero with 100% errors 
and does not show clear signal of a linear dependence on t. We extract the tensor charge 
both for the connected and disconnected contributions fitting the data of Rrqif) with a 
linear form (iTf) over 6 < t < 11. 



(20) 



The results for 5q are corrected by the one- loop renormalization factor |T6| , pO 

ZT,{^ia,a^{l/a)) = [l - ^) [l - (^|- ln(H + 0.44^ ajj^{l/a) 

The tensor charges 5u, 6d and 6s in the MS scheme are tabulated in Table ^. 

In Fig. 1^ we plot the tensor charges as a function of the lattice bare quark mass rriud- 
The connected contributions for u and d quarks decreases in magnitude as the quark 
mass decreases. This trend is the same as for the axial charge. For the disconnected 
contribution our result is consistent with zero, which indicates that the OZI violating 
effects are negligible for the tensor charge. 
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Figure 4. Same as Fig. ^ for Rx^it). Solid Figure 5. Same as Fig. |^ for the tensor 
lines are linear fits over 6 < t < 11. charge. 



To extract the physical tensor charge for u and d quarks the fitted values are linearly 
extrapolated to the chiral limit iriud = 0. For the strange quark contribution we first make 
linear interpolation to the physical strange quark mass m^a at each fixed value of K^d 
for nucleon, and then the results are linearly extrapolated to the chiral limit. The final 
results for the tensor charge is 

5i: = 5u + 5d + 5s = +0.839(60) - 0.231(55) - 0.046(34) = +0.562(88), (21) 

where the error of of the sum of u, d, s contributions are estimated by quadrature. 
Due to the smallness of the disconnected contributions, the flavor singlet tensor charge 
is not much suppressed from its quark model value. This is in contrast to the flavor 
singlet axial charge AE which suffers from large suppression. The ratio |5g/Ag| extrap- 
olated to the chiral limit gives |5u/Am| = 1.32(15) and \5d/^d\ = 0.67(18), which are 
quantitative different from the prediction of the nonrelativistic quark model \6q/Aq\ = 1 
{q = u, d) and also have qualitative difference from the prediction of relativistic quark 
models \6q/Aq\ > 1 {q = u,d) mentioned before. The smallness of Sqdisc could be related 



Table 3 

Proton tensor charges at /i = 1/a in the MS scheme as a function of Kud- Sd^isc equals 





Values at Kc 


are obtained by a linear fit 


in 1/K. 
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'^^conn 
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ir,=0.1665 


0.1640 


0.1600 


0.1600 


1.072(50) 


-0.251(21) 


-0.008(12) 


-0.022(16) 


-0.013(14) 


-0.008(12) 


0.1640 


0.994(11) 


-0.220(52) 


-0.027(28) 


-0.041(33) 


-0.027(28) 


-0.013(25) 


0.1665 


0.948(31) 


-0.196(15) 


-0.044(69) 


-0.044(69) 


-0.035(58) 


-0.009(51) 


Kc 


0.893(22) 


-0.180(10) 


-0.054(54) 
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to the C (charge conjugation)-odd and chiral-odd nature of the tensor operator qa^^l^fl 

Possible systematic errors originates from the scahng violation effects and the quenched 
approximation. We may estimate that the magnitude of these systematic errors is ~ 
20% as for the axial charges. Toward a definitive determination of the tensor charge a 
repetition of the calculation including dynamical quark effects with smaller lattice spacing 
is required, which we leave for future investigations. 
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